Nutation versus angular dependent NQR spectroscopy and the impact of underdoping 

on charge inhomogeneities in YBa2Cu 3 O y 
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We describe two different nuclear quadrupole resonance (NQR) based techniques, designed to 
measure the local asymmetry of the internal electric field gradient rj, and the tilt angle a of the 
main NQR principal axis z from the crystallographic axis c. These techniques use the dependence of 
the NQR signal on the duration of the radio frequency (rf) pulse and on the direction of the rf field 
Hi with respect to the crystal axis. The techniques are applied to oriented powder of YBa2Cu30 y 
fully enriched with 63 Cu. Measurements were performed at different frequencies, corresponding 
to different in-plane copper sites with respect to the dopant. Combining the results from both 
techniques, we conclude that oxygen deficiency in the chain layer lead to a rotation of the NQR 
main principal axis at the nearby Cu on the Cu02 planes by a ~ 20° ± 5°. This occurs with no 
change to n. The axis rotation associated with oxygen deficiency means that there must be electric 
field inhomogeneities in the CUO2 planes only in the vicinity of the missing oxygen. 
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I. INTRODUCTION 

Hole doping the CuC>2 planes of the cuprates creates 
natural inhomogeneities. Determining the origin and 
form of these charge inhomogeneities is essential for our 
understanding of their transport properties. There are 
several theoretical works suggesting that phase separa- 
tion is a characteristic of the Hubbard model, and there- 
fore an intrinsic property of the CuC>2 planes. 1 Moreover, 
it has been suggested that the phase separation, possibly 
in the form of strips, can explain many of the unusual 
properties of the cuprates and leads to the superconduc- 
tivity (for example, Ref.— ). 

Indeed, some of these materials show phase separation 
with segregated hole-rich and hole-poor regions. Such in- 
homogeneities were found in I^-ajSr^CuiC^ by several 
methods such as nuclear quadrupole resonance (NQR), 3 
neutron scattering^ and muon spin relaxation (/iSR). 5 
In addition, Scanning tunnel microscope (STM) experi- 
ments on underdoped Bi2Sr2CaCu20j, (Bi-2212) showed 
local density-of-states modulations^ and inhomogeneity 
of the superconducting gap on the samples' surface that 
can be associated with the distribution of holes in the 
planes X 

In YBa2Cu30 y (YBCO), phase separation was found 
at very low doping levels, up to YBCO6.35, with 
neutron scattering from phonons related to charge 
inhomogeneity^ /iSR measurements show the existence 
of a spin-glass phase for a similar doping ranged The 
highest doping in which magnetic order was found in 
YBCO was in YBCOg.e 10 using neutron scattering. 89 Y 
NMR study in YBCO a for y = 7 and 6.6 showed no 
phase separation at allM Direct detection of charge in- 
homogeneity in YBCOj, is nearly impossible since STM 
measurements are very difficult due to oxygen loss in vac- 
uum, and surface problems. Various STM experiments in 
this compound showed an inhomogeneousi^ or relatively 
homogenous surface,— depending on the surface prepa- 
ration procedure. In these circumstances NQR might be 
the only tool to probe charge inhomogeneities in YBCO. 



Moreover, NQR experiments are sensitive to the charge 
distribution in the bulk and not just on the surface. 

In this paper, we present NQR measurements on 
63 Cu enriched YBCO. This YBCO compound has nar- 
row NQR resonance lines, which allows us to distin- 
guish between different in-plane copper [Cu(2)j and asso- 
ciate each line with a local environment. We then apply 
two additional NQR based methods called nutation spec- 
troscopy and angle dependent NQR (ADNQR). The two 
methods are set to measure the quadrupole interaction 
asymmetry parameter r\ and the tilt angle a of the NQR 
main principle axis z from the crystal axis c. This allows 
us to measure charge inhomogeneity in the Cu02 plane. 
The main advantage of these techniques is that they can 
be employed at every frequency in the NQR spectrum. 
Therefore, unlike NMR, rj and a can be determined for 
each ionic environment separately. Our main conclusion 
is that for Cu(2) in unit cells with no oxygen deficiency 
a ~ and r\ ~ 0. However, for every Cu(2) in a unit cell 
where an oxygen is missing, and only in these unit cells, 
a ~ 20° and r\ ~ 0. Such a tilt of z from the crystal axis 
c in special unit cells implies charge inhomogeneities in 
the Cu02 plane only next to oxygen deficiency. 

The paper is organized as follows. In Sec. [IT] we give a 
general description of NQR and then introduce the less 
known experimental methods: the nutation spectroscopy 
and the ADNQR. Theoretical calculations for these tech- 
niques are also shown. In Sec. IIHI we describe the exper- 
imental setup. The experimental results are presented 
in Sec. IIVI Finally the summary and conclusions on the 
experimental techniques and their application to YBCO 
are presented in Sec. El 



II. THE EXPERIMENTAL TECHNIQUES 
A. NQR 

In an NQR experiment one uses an rf magnetic pulse in 
order to cause transitions between the nuclear energy lev- 
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els. Nuclei with spin I > 1/2 can be viewed as positively 
charged oval objects. As a result, their energy inside 
a solid depends on their orientation in the electrostatic 
potential V(r) generated by the other nuclear and elec- 
tronic charges. When the nuclear poles are close to posi- 
tive charges their energy is high, and when the poles are 
close to negative charges the energy is low. The energy 
difference between different orientations is determined by 
the electric field gradient (EFG) tensor Vij — d ®.g x a * 
the position of the Cu nuclei. The directions can be cho- 
sen so that Vij is diagonal. These directions are known as 
the principal axis of the EFG. Due to Laplace equation 
(V X x + V yy + V zz — 0) the NQR Hamiltonian is deter- 
mined by only two parameters, v q and 77, and is given 
by: 



I 2 +V(I 2 X -Iy)} 



(1) 



where I a are the nucleus spin operators; v q is a frequency 
scale, determined by the main EFG component V zz , and 
the quadruple moment of the nucleus Q; 
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(2) 



is the asymmetry parameter with the convention | V zz | > 
\Vyy\ > \Vxx\, and therefore < 77 < 1. 

For the spin 3/2 copper nucleus, the NQR Hamiltonian 
has two doubly degenerate energy levels, and therefore 
only one resonance frequency, given by: 



NQR 



1+ V 



(3) 



Since v q is proportional to V zz , it holds information on 
the charges surrounding the nucleus. Therefore nuclei 
sitting in different electronic environments will have dif- 
ferent resonance frequencies. 77 is a parameter that holds 
information on the anisotropy of the surrounding charges 
in the x-y plane. In the case of axial symmetry, 77 = 0. 
It is clear from Eq. that in the standard NQR exper- 
iment where only the resonance frequency is measured, 
77 and v q cannot be separately determined. However, the 
properties of the rf field turn out to be handy. 
The rf pulse Hamiltonian is 

Hrf = r yhH 1 ■ Isinwt, (4) 

where, in general, the rf field can be described as 



Hi = H t T 



[sin(0) cos(</5), sin(0) sm(<j>), cos(0)] 



(5) 

where 9 and cf> the polar angles relating the coil axis to 
the quadrupolar reference frame. 

The signal in the coil after a time t from an 011- 
resonance rf pulse (uj — ujq = 2n Jnqr) with a duration 
t p is given bj*^ 



I(t p , t, 9, <j>, 77) cx A(6>, (f>, 77) sin[A(6', 0, rj)u)ib p ] sin(w Q i) 



(6) 



where u)\ — jHi (7 is the nucleus gyromagnetic ratio), 
and A is an angular factor given by 



A(0, <p, 77) = ^ rial + rial + r 2 z a? z , 



where 



T) 



= (77 + 3,77-3,277)- 



(7) 



(8) 



Hence, the signal at time t is a function of both the du- 
ration of the rf pulse and the polar angles between the 
rf direction and the EFG. In the following sections we 
review two methods that use these parameters to deter- 
mine 77, and the angle a between z and c. 



B. Nutation spectroscopy 

The nutation spectroscopy NQR, developed by Harbi- 
son et aZ . 15 i 16 is a method in which the signals are mea- 
sured as a function of the duration of the rf excitation 
pulse t p . This method is used to determine the asymme- 
try parameter 77. Integrating Eq. ([6]) over the angles for 
isotropic powder and Fourier transforming with respect 
to t p gives a powder pattern line shape that is described 
by 

t, 77) oc / dej) I sm9d9\(0,(j),r)) 



x / e l " p p dt p sin(A(0, 0,77)^1 tp)sin(wgt). 



(9) 



The Fourier transform part of this integral, over t p , gives 
a delta function at the angular frequency, 



u p (9,<t>,ri) = \(9,<P,t])uj 1 



(10) 



The angular integrals required to complete the calcula- 
tion of I{u)p,t,rj) demand numerical integration. Nev- 
ertheless, the peaks of this function can be obtained by 
looking at the extremum of u> p (9, 4>, rf) with respect to 9 
and <j). This gives three singularities 
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In the special case of 77 = we find lujj — ujjjj and 
u>i = 0, and therefore there is only one sharp frequency 
at 0.866^1. As 7/ grows the three frequencies separate into 
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FIG. 1: (Color online) Calculated on-resonance nutation spec- 
tra for: (a) powder samples, (b) Perfectly oriented powder 
where z||c. (c) Oriented powder where z has up to 20° an- 
gle from c (d) Perfectly oriented powder where z||b. For the 
oriented powders the rf transmission is perpendicular to the 
c direction. 
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FIG. 2: (Color online) The unit cell of YBCO. An illustration 
of the rotation of the EFG main principal axis from the crys- 
tallographic directions, due to the missing oxygen atoms in 
the chains. The angle between z and c is a. The x direction 
is in the a-b plane making an angle ip with a 



the three singularities and the pattern broadens. From 
these singularities r\ can be extracted by r\ — (ujuj — 
ujii)/{ujiii+loii). FigureQTa) shows the powder patterns 
for different values of r\ for randomly oriented powder. 

In the case of an oriented powder, the crystallographic 
c direction and the EFG z direction are parallel but the 
crystallographic a (and b) directions are random with 
respect to the EFG x (and y) directions. The rf trans- 
mission is done perpendicular to c (and z). Figure [Tfb) 
shows the nutation line shapes in this case. There is no 
big difference between powder and oriented powder for 
rj = but as r\ grows the difference between the two 
cases becomes clear. 

In general, however, z does not have to be exactly par- 
allel to the c direction of the lattice. This situation is 
demonstrated in Fig. [5] The angle between z and c is 
a. The x direction is assumed to be in the a-b plane 
making an angle if with a. Figure [UJc) shows the line 
shapes when a = 20° and ip is averaged. One can see 
that with the nutation spectroscopy it is very difficult to 
distinguish between the case of perfect orientation (panel 
b) and partial orientation (panel c) . Finally, in Fig. [IJd) 
we present a case of perfectly oriented powder but with 
z parallel to b. This case is fundamentally different from 
previous cases especially when 77 = 1. 

For a spin-echo sequence, after a time r from the first 
pulse there is a refocusing pulse with duration t r . This 
adds an additional factor of sin 2 [A(6*, <fi, r))ujit r ] to the 
signal. In an NQR experiment, since the spin rotation 
frequency depends on the orientation of the lattice with 
respect to the coil, the second pulse cannot perfectly re- 
focus all the magnetization. Harbison et a/. 16 showed 
that the additional factor does not change the nutation 
frequencies ujijijii; however, it does change the relative 
intensities of these frequencies with respect to each other. 

The main advantage of the nutation method is that it 
is relatively simple to execute; it can be carried out on 
a simple NQR spectrometer without a static magnetic 
field or additional modifications. This method can also 
be implemented for powders. It allows the determination 
of 77 at every point of the NQR spectrum (unlike NMR, 
where 77 can be determined only from the entire spectrum 
with no local resolution). 

C. Angle dependent NQR 

The angle-dependent NQR technique on an oriented 
powder was originally developed by Levy and KerenJi In 
this technique the signal intensity for a given frequency 
is measured as a function of 9, the angle between the di- 
rection of the rf field and the crystallographic c direction. 
This technique is sufficient to determine 77 when the EFG 
z direction is in the crystallographic c direction. As we 
mentioned before this does not have to be the case. In the 
ADNQR experiment the sample is rotated with respect 
to the symmetry axis of the coil (see inset of Fig. [3]) , and 
the echo intensity is registered for each angle. 
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FIG. 3: (Color online) Theoretical echo intensity curves for 
various values of 77 calculated from Eq. (115[l , 8 is the polar 
angle between the rf field and the c direction of the lattice, in 
an oriented powder. In this case the z direction of the EFG 
is assumed to be parallel to c. Inset: basic angle-dependent 
NQR configuration. A sample with a preferred direction is 
inserted into the coil. The angle 6 between them can be varied 
with a motor. 
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FIG. 4: (Color online) Theoretical echo intensity curves for 
77 = and various values of a calculated from Eq. (|19[) . 9 is 
the polar angle between the rf field and the c direction of the 
lattice, in an oriented powder. 



From Eqs. {TJ), ([3]), and one can see that for 77 = 0, 
H.q, and TL r f commute when 8 = 0. In this case there 
will be no spin transitions and no signal. For 77 > 0, H.Q, 
and Ti r f do not commute even for 8 — 0. In this case 
we expect a signal even when Hi is in the z direction. 
In general, Levy and Kerenii showed that for a tt/2 — 
t — 7r pulse sequence, when Hi is given by Eq.( [5]), the 
magnetization in the coil at the time of the echo is given 
by 



\(d,<j),r})tikjQ 
2K B T 



sin 3 (\{d,(j>, 77)^/2), (14) 



where t^n is the duration of the n/2 pulse, and A(#, 0, 77) 
is the same as in Eq. ([7]). 

In the case of an oriented powder, again with z and 
c parallel, and the a (and b) random, M is obtained by 
averaging over <j), namely, 



M{n,6) = ±- [ M{o,<p ir] )d4> 

27r Jo 



(15) 



Theoretical echo intensity curves as a function of 8 for 
various values of 77 are presented in Fig. [3] A fit of ex- 
perimental data to these theoretical curves can give the 
value of 77. 

Next we discuss the case where z has an angle a with 
c as in Fig. [2] In this case TLq and Ti. r f again do not 
commute even for 8 = and we expect a signal even 
when Hi is in the z direction, for all values of 77. To 
evaluate the signal intensity we can take the unit vector 
pointing in the direction of Hi in the crystal coordinate 
system to be r = xsin# + zcos0 and express it as ?■ in 
the EFG coordinate system using the Euler angles^. If 



we assume, again for simplicity, that x is in the a-b plane 
with an angle tp from a as in Fig. [5] then, 



r {a, ip, 



and 



1 

cos a sin a 
— sin a cos a 



Ha, Lp, 8,i]) = 



cos ip sin ip "1 
— sin ip cos f r (8) 
1 

(16) 



(17) 



\J ri {a, ip, 8)a 2 x (77) + rf {a, ip, 8)a 2 y (77) + fi 2 (a, ip, 8)a? z (77) . 
In the case of 77 = 0, 

\(a,ip,9,r) = 0)= (18) 

V3 



y cos 2 ip sin 2 8 + (sin a cos 8 — cos a sin ip sin 8) 2 

and the magnetization in the coil as a function of the 
rotation angle 8 is now 



M(a,6,T] = 0) = 



2k B T 



(19) 



2tt 



x 7j- / A(a, ip,8,r\ = 0) sin 3 [X(a, ip, 8, 77 = 0)w 1 i w / 2 ] dip. 



Figure 2] shows the magnetization in the coil according 
to Eq. (jTUJ) for different values of alpha. 
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III. THE EXPERIMENTAL SETUP 

The NQR spectra were obtained using a spectrome- 
ter with a home-made automated frequency sweep. The 
measurements were performed in a coil tunable from 25 
to 33MHz. The data were obtained by applying a spin- 
echo sequence. Both for the nutation and the ADNQR 
techniques it is of great importance for the rf field to be 
homogeneous. For the nutation experiment we measured 
samples with a small volume inside a long cylindrical coil. 
The coil's wire goes through a current monitor. This cur- 
rent monitor allows us to perform measurements at dif- 
ferent frequencies with the same current through the coil 
and therefore the same Hi. For the ADNQR we used 
a spherical coil that gives a more homogenous magnetic 
field with a better filling factor. The experiment is fully 
automated, the sample holder is connected to a motor 
which rotates the sample and can be controlled from the 
computer. 

The measurements were done on YBCOj, oriented pow- 
ders, with different doping levels. The data were taken 
at 100 K, which is above T c for all samples. Previous 
Cu NQR measurements on YBCO (see RefJ£) were per- 
formed on samples containing both Cu isotopes, Cu 63 
and Cu 65 , so the frequency lines consisted of doublets of 
Cu 63 -Cu 65 . In this work, for a clearer understanding of 
the NQR signals, the YBCO samples were all made of 
enriched copper, meaning that these samples contained 
only the Cu 63 isotope. This allows us to distinguish be- 
tween the different contributions to the NQR line from 
different local environments. 



IV. RESULTS 
A. Cu NQR lines of enriched YBCO samples 

The frequency sweep lines of YBCO y samples with dif- 
ferent doping at 100 K are presented in Fig. [5] The spec- 
trum is normalized by f 2 in order to correct for popula- 
tion difference and the induced signal in the coil. For the 
fully doped sample only one resonance is observed. Be- 
tween y — 6.85 and 6.73 two resonances are observed. At 
y = 6.68 there are clearly three resonances. The middle 
one disappears upon further oxygen reduction and below 
y = 6.45 again only two resonances are seen. For each 
sample, the resonance frequencies were extracted from 
the spectrum. In Fig. Ufa) we plot the oxygen doping as 
a function of these frequencies. From this plot it is clear 
that there are three different Cu sites. 

Earlier NQR measurements done on YBCO y with both 
Cu isotopes present, found similar resonance frequencies 
for similar y values . 19 i 20 We therefore assign the different 
resonance frequencies of each sample to different ionic 
environment, based on these earlier work. For high dop- 
ing level (y > 6.5) all signals come from the plane Cu(2). 
There are three different types of environment that affect 
the Cu(2) resonance frequency. They were classified in 
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FIG. 5: (Color online) NQR frequency sweep on YBCOy. The 
solid lines are Gaussian fits performed in order to determine 
the resonance frequencies. 



terms of the number of oxygen surrounding the chain cop- 
per Cu(l) neighboring the detected Cu(2). Cu(l)4 stands 
for a Cu(2) next to a full chain as in YBCO7; in this case 
the frequency is the highest. Cu(l)3 means Cu(2) whose 
neighboring Cu(l) is missing one oxygen; the middle fre- 
quency belongs to this Cu(2). Finally, Cu(l)2 is when 
the neighboring chain to the Cu(2) is empty; the lowest 
frequency belongs to this case. The three possible envi- 
ronments of the Cu(l) are shown schematically in Fig. [7] 

It is clear from Fig.[5]that as the doping decreases, the 
amplitude of the lower frequencies resonances increases at 
the expense of the high-frequency resonance. This makes 
sense with our line assignment, as at lower doping there 
are less Cu(2)4 sites and more Cu(2) 3 and Cu(2) 2 sites. 
In addition, at oxygen doping close to 6.5 (i.e., 6.56 and 
6.45) the signal from Cu(2)3 is very hard to detect. This 
is due to the formation of the ortho II phase at these 
dopingsi^i 

In Fig. [6ja) we depict the oxygen level y vs. Jnqr- 
Clearly, the more holes are introduced into the planes the 
more shifted is the frequency. In Fig. [Hfb) we plot the 
superconducting transition temperature Tc as a function 
of Jnqr- For each site, the lower Tc the more shifted 
are the peaks toward low frequencies. It seems that the 
rate of change of Tc with Jnqr] dTc / df^QR is identical 
for all sites as demonstrated by the solid lines. Assum- 
ing that the NQR resonance frequency is proportional to 
the number of charges (which holds in the case where r\ 
is constant), this implies that Tc is proportional to the 
number of free charges in the plane. Since Tc is also pro- 
portional to the superfhrid density n s ^ this means that 
the rate of conversion of holes to superconducting holes 
is constant. 
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FIG. 6: (Color online) (a) The oxygen doping for different 
YBCO samples as a function of the NQR frequencies, (b) The 
superconducting transition temperature of these samples as a 
function of the NQR frequencies. The blue circular symbols 
are resonance frequencies coming from the Cu(l) site. The 
red lines are a fit to three parallel lines, for the three different 
Cu(2) sites. 
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FIG. 7: (Color online) Schematic illustration of the Cu site 
in YBCO with locally different oxygen coordinations. 



B. Nutation spectroscopy 



In the nutation experiment we used a pulse sequence 
shown in the inset of Fig. [5J The measurement is done 
by applying this sequence at the resonance frequency of 
each site, with varying excitation pulse length t p . The 
entire echo is integrated and recorded as a function of 
t p . An example for raw data for such experiment on a 
sample with y=7-S is shown in Fig. [8ja). The Fourier 
transform over t p of this data is presented in Fig.[8^b). 

For the Cu(2) of fully doped YBCO it was 
established 2 ^ 2 ^, that the principal component of the EFG 
is in the c direction of the lattice and that 77 = 0. Looking 
at Eqs. (|10jl - (fT3|) . in the case of i] = there is only one 
nutation frequency and the ratio tOp/ui is 0.866. There- 
fore, from the nutation frequency of the Cu(2) resonance 
of YBC07_5 we can extract uj\. Since we worked with 
a constant rf field H±, we normalized the frequency axis 
for all samples by w\. We then apply this technique to 
three YBCOj, samples. The results are shown on the left 
panels of Fig. [9] For each sample the measurements were 
done at the frequencies marked with arrows on the NQR 
spectrum, extracted for clarity from Fig. [5j and shown in 
the right panels. 

Figure [9]^a) shows the nutation spectrum for 
YBC0 7 _,5, measured in both the Cu(2) (31.5MHz) and 
the Cu(l) (22MHz) resonance frequencies. After nor- 
malizing the frequency axis by w\, as explained above, 
the Cu(2) has a sharp nutation frequency with u> p /u>i = 
0.866. However, Cu(l) has a much broader spectrum 
with uj p jjj/uji = 0.52 and u p u which is difficult to deter- 
mine. The result for the Cu(2) is consistent with 77 = 
[see theoretical nutation spectra in Fig. [TJb)] . The re- 
sult for Cu(l) is consistent with r\ = 0.95 ± 0.05, and 
z pointing in the b direction, as shown in Fig. QJd). 
These results are in agreement with NMR results on 
YBC0 7 _ 5 (Refs322) that measured ?7~0 for Cu(2), 
and 77 ~ 1 for Cu(l). Similar nutation experiments on 
powder YBC07-5 at room temperature were done by 
Vega. 24 His results are 77 = for the Cu(2) and 77 = 0.8 
for the Cu(l). 

Figures. Mb) and^c) show the nutation spectroscopy 
results for YBCO y with lower doping levels. For these 
samples we measured only at the resonance frequencies of 
the Cu(2). The sample with y = 6.73 has two resonance 
lines for two different Cu(2) ionic environments, and the 
sample with y = 6.68 has three Cu(2) environments. The 
nutation experiment shows that for these samples, for 
all three different types of Cu(2) ionic environments, the 
nutation spectrum is not different from that of the y—7- 
S case, therefore 77 ~ 0. This conclusion was verified 
by using several different H\ and refocusing pulse time 
t r . However, as mentioned in section Hi B[ the nutation 
experiment is not sensitive enough to small differences in 
the orientation of the EFG [see Figs.QJb) andQJc)]. 
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FIG. 8: Inset: the pulse sequence of a nutation experiment, 
(a) An example of raw data of a nutation experiment for y=7- 
5 at 31.5MHz. The refocusing pulse t r was 4.4/is and r was 
32/^s. (b) The data after Fourier transform. 
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FIG. 9: (Color online) (a) Nutation spectra for YBC07_i 
with natural abundance of 65 Cu and 63 Cu. [(b) and (c)] Nuta- 
tion spectra for YBCO6.73 and YBCO6.68 enriched with 63 Cu. 
[(d)- (f)] Cu NQR line shape for these three samples, extracted 
from Fig. [5] The arrows show the frequencies where nutation 
spectroscopy is applied. 



C. ADNQR 

As mentioned before, in a fully doped YBCO z is par- 
allel to the c axis of the lattice. Therefore the ADNQR 
technique can be applied to this sample, using Eq. (fT5)) . 
The result indeed gives 77 — as was published in Ref. 25 
and in agreement with the nutation experiment described 
above. For lower doping levels however, the exact direc- 
tion of z is unknown. In a previous publication in Ref. 25 
we also performed an ADNQR experiment on a YBCO 



FIG. 10: (Color online) Right panels: Cu NQR line shape for 
YBC07-* with natural abundance of 65 Cu and 63 Cu, and for 
YBCO6.68 enriched with 63 Cu. The arrows show the frequen- 
cies where ADNQR is applied. Left panels: the echo intensity 
as a function of the angle 9 between the rf field and the c di- 
rection of the lattice, for the two samples. The solid lines are 
fits to Eq. (PH . 



sample with doping of y = 6.68. Due to lack of infor- 
mation we assumed that zj|c at all doping. The experi- 
mental results from Ref.— are shown by the dotted line 
in Fig. [TW b). Fitting these data to Eq. ([15)) gave a high 
value of 77 in contradiction to the nutation experiment 
findings in Fig. [9^c) of the present work. The way to 
settle this contradiction is by allowing z to be tilted with 
an angle a from c as in Fig. [5] As mentioned above, in 
this case the NQR signal does not disappear at 8 = 0, 
even in the case of rj = 0, as result of z EFG component 
in the a — b plane. 

Figures. [TI|a) and [TUKb) show the ADNQR data for 
the two YBCO samples with doping of y—7-S and 6.68. 
The solid lines are new fits, this time to Eq. (fT9|) . namely, 
i] = and tilted principle axis. The fit allows for a fi- 
nite base line for each sample to account for some un- 
known amount of misalignment. This misalignment is a 
result of the non-perfect orientation and some inhomo- 
geneity of the induced field in the coil. These samples' 
frequency spectrum, taken from Fig. O are presented 
again in Figs. ITOT c) and [TOM) for clarity. The arrows 
in this figure mark the frequencies where ADNQR was 
applied. 

For the main peak of y=7-S the best fit is obtained 
with a = and rj = 0, as expected. The same result was 
obtained at the high frequency of the y = 6.68 sample, 
which is also associated with Cu(2) in an environment 
of full chains; Cu(l)zi. In contrast, for the other oxygen 
environments the signal at 9 — > is clearly above the 
background. Keeping in mined our nutation spectroscopy 
results, showing that 77 — 0, this suggests that the value 
of a is larger than zero. The fit to Eq. (TIT?)) gives a ~ 
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20° ± 5°. This is the main finding of this work. 

V. SUMMARY AND CONCLUSIONS 

We used nutation spectroscopy and ADNQR to mea- 
sure the quadrupole interaction asymmetry parameter 
r\ and tilting angle a of the main component of the 
EFG principal axis z from the crystallographic c direc- 
tion. Calculations show that nutation spectroscopy pow- 
der patterns for oriented powder are very sensitive to the 
value of r], but they are not very sensitive to a for a up 
to 20°. In contrast, the ADNQR line shape is sensitive 
to both r\ and a. In other words, ADNQR is sensitive 
to any breaking of the axial symmetry. However, it is 
difficult to determine both a and 77 from ADNQR alone. 
The combination of the two methods can give a very good 
estimate on both parameters. We implement these tech- 
niques on the in-plane Cu atom of the YBCO compound 
with different doping levels. We conclude that for both 
YBC07_i and YBCOg.68, V is approximately zero for 
all Cu(2) resonance frequencies. However, the ADNQR 
showed that for the underdoped sample there is a clear 
difference between the Cu(2) neighboring a Cu(l) in a 
full chain as in a fully doped sample [Cu(l)4] and the 
Cu(2) neighboring a Cu(l) in an empty or a half filled 
chain [Cu(l)2-Cu(l)3]. This difference can result from a 
small change in the EFG principal axis z with respect to 
the lattice directions c. 

The motivation for these experiments was to measure 
possible charge inhomogeneity or electronic phase sepa- 
ration in the YBCO compound. Both the nutation spec- 
troscopy and the ADNQR for the Cu(2) from a fully 
oxygenized local environment Cu(l)4, even for lower av- 
eraged doping, show a homogeneous charge distribution 
in the plane. Our combined experiments also imply that 
only for Cu(2) neighboring a Cu chain with missing oxy- 



gen Cu(l)2 or Cu(l) 3 is the principle axis z of the EFG 
tensor not along c, namely, a ~ 20° ±5°. Therefore, only 
from the point of view of a Cu(2) in YBCO6.68 neighbor- 
ing a Cu(l)2 or Cu(l)3, are the a and b crystallographic 
directions less identical than in YBCO7. This means that 
any charge inhomogeneity in the plane is correlated di- 
rectly with the O dopant atoms, and therefore cannot 
be an intrinsic property of the Cu02 planes alone. It is 
either generated or pined by missing oxygen. 

McElroy et al^ came to a similar conclusion by per- 
forming spectroscopic imaging scanning tunneling mi- 
croscopy on Bi-2212 samples. They found strong cor- 
relation between the position of localized resonance at 
-960meV identified with interstitial oxygen dopants and 
the size of the local spectral gap. To understand these re- 
sult, Nunner et al£l presented a theoretical model where 
the dopants modulate the pair interaction locally on an 
atomic scale. They calculated the correlation between 
the local density of states and the dopant modulated 
pair interaction potential. They showed that this model 
agrees with McElroy's experimental results on Bi-2212. 
A more recent theoretical work by Mori et al£2- identified 
two mechanisms by which the position of the apical oxy- 
gen atoms can modulate the pairing interaction within 
the Cu02 planes. 

Our result for the YBCO compound reinforces the sur- 
face experiments done on Bi-2212. It shows that the cor- 
relation between the electronic spatial variation in the 
plane and the dopant exists not only in Bi-2212 and it is 
a property of the bulk and not only of the surface. 
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